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(54) PICTURE PROCESSOR 
(57)Abstract: 

PROBLEM TO BE SOLVED: To provide a picture processor for operating quick 
distortion correction. 

SOLUTION: A distortion correcting part 203 of a picture processor 202 generates a 
picture after correction from a picture before correction including geometrical 
distortion generated due to the characteristics of an X-ray photographic system for 
correcting distortion generated due to the characteristics of an image intensifier. 
Therefore, a picture after correction can be generated by the following procedures, 
that is, (1) the measurement of rough distortion distribution (measuring means), (2) 
the calculation of fine distortion distribution (calculating means), and (3) the decision 
of a concentration value (concentration value deciding means). 
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CLAIMS 



[Claim(s)] 

[Claim 1] In the image processing system which creates the image after amendment 
from the image before amendment including a geometric distortion which photoed the 
phantom which has specific geometric structure by the photography system, 
originated in the property of the photography system concerned, and was produced A 
measurement means to measure the coarse strain distribution which shows the 
distorted point of the specifying point specified according to the geometric structure 
of said phantom from the photography image of said phantom, A calculation means to 
compute the fine strain distribution which divides the image after said amendment into 
two or more minute fields, and shows the distorted point of the representation point in 
the minute field concerned using said coarse strain distribution, The image processing 
system characterized by providing a concentration value decision means to determine 



the concentration value of the minute field on the image after said amendment based 
on the concentration value of the minute field on the image before said amendment 
specified by the distortion place of said representation point. 

[Claim 2] Said calculation means is an image processing system according to claim 1 
characterized by computing the fine strain distribution which shows the distorted 
point of the representation point concerned by carrying out the 2-dimensional 
interpolation operation of the coarse strain distribution which shows the distorted 
point of at least four specifying points near [ in said minute field ] the representation 
point. 

[Claim 3] Said calculation means is an image processing system according to claim 1 
characterized by providing a means to divide the image after said amendment into two 
or more minute fields to which the core of an image differs in magnitude from a 
periphery, and a means to compute the fine strain distribution which shows the 
distorted point of the inner each representation point of said minute field based on 
said coarse strain distribution. 

[Claim 4] A setting means to set up an area of interest on the image after said 
amendment is provided further. Said calculation means A means to divide the image 
after said amendment into two or more minute fields from which magnitude differs the 
inside of the field equivalent to the area of interest set up by said setting means, and 
out of the field, The image processing system according to claim 1 characterized by 
providing a means to compute the fine strain distribution which shows the distorted 
point of the inner each representation point of said minute field, based on said coarse 
strain distribution. 

[Claim 5] A means to pinpoint the minute field on the image before said amendment 
equivalent to the minute field on the image after said amendment according to the 
distortion place of the representation point that said concentration value decision 
means was computed by said calculation means, The image processing system 
according to claim 1 characterized by providing a means to determine the 
concentration value in the minute field on the image after said amendment, based on 
the concentration value in the minute field on the image before said amendment 
specified by this means. 

[Claim 6] It is the image processing system according to claim 1 which said phantom 
has grid-like structure, and said measurement means integer-izes the distance during 
the lattice point of said phantom on the image after said amendment, measures the 
coarse strain distribution which shows the distorted point of each lattice point when 
the integer-ization concerned was made, and is characterize by for said calculation 
means to divide the image after said amendment so that it may be locate at the core 
whose representation point of two or more of said minute fields is a pixel. 
[Claim 7] The image processing system according to claim 6 characterized by 
providing further a means to expand or reduce the image after said amendment by 



changing the distance during said lattice point. 

[Claim 8] The image processing system according to claim 7 characterized by 
distance between said grids into the nearest integer. 



DETAILED DESCRIPTION 



[Detailed Description of the Invention] 
[0001] 

[Field of the Invention] This invention relates to the suitable image processing system 
for distorted amendment of the image outputted from various modality, such as an 
image processing system which amends distortion of an image, X-ray picture 
photography equipment which uses an image intensifier (I.I.) especially, or a 
computerized-tomography scanning equipment. 
[0002] 

[Description of the Prior Art] While this detects the transparency X-ray of analyte in 
the X-ray-diagnosis equipment used widely in current and a medical site, using an 
image intensifier as an X-ray detection system, what changes into an optical image 
and photos the obtained optical image by the TV camera system is known. There are 
digital radiography and fluorography by X-ray image intensifier as the roentgenography 
approach of digitizing the picture signal collected with this equipment. The video signal 
from a TV camera is changed into a digital signal by the analog / digital (A/D) 
transducer, and in order to perform various image processings, it is inputted into a 
computer. 

[0003] The image intensifier has the structure which piled up an X-ray fluorescence 
side and the photoelectric surface. The X-ray which penetrated analyte is inputted 
into an X-ray fluorescence side, when the light from this X-ray fluorescence side 
stimulates the photoelectric surface, a photoelectron is emitted and image formation 
of the visible image (optical image) is carried out to an output phosphor screen. 
[0004] There is a problem that the geometric distortion of spool type distortion, S 
character-like distortion, etc. will be included in the output image in such an image 
intensifier. For example, photography of an object object (phantom) as shows a wire to 
drawing 12 (a) arranged in the shape of a tetragonal lattice on the flat surface will 
output the image distorted as shown in drawing 12 (b). 

[0005] For example, when performing the image processing for extracting useful 
information from an output image to a diagnosis in the above-mentioned digital 
radiography, such distortion can become a big factor with error. In image processings, 
such as clinical analyses, such as rate count of the vasoconstriction, and 
ventriclus-sinister volume count, or three-dimensions reconstruction, the effect of 



the error which originates distorted becomes remarkable especially. 

[0006] In view of such a background, also in order to perform image processings, such 

as clinical analysis (bibliography ''the Toshiba digital 

fluorography-by-X-ray-image-intensifier system CLINICAL ANALYSIS", Toshiba 

Corp., and Toshiba Medical Systems Co., Ltd.") or three-dimensions reconstruction, 

appropriately and with high precision, distorted amendment of the image resulting 

from an image intensifier is indispensable, and distorted amendment is performed in 

the phase of pretreatment in these image processings by the former. 

[0007] However, in the distorted amendment currently performed from the former, 

since fine strain distribution is calculated for every pixel and he is trying to obtain an 

amendment image, there is a trouble that processing takes immense time amount. At 

a clinical place, a highly precise diagnostic image is offered quickly, and also in order 

to aim at improvement in diagnostic ability, it is anxious for improvement in the speed 

of the distorted amendment in an image processing. 

[0008] 

[Problem(s) to be Solved by the Invention] Therefore, the purpose of this invention is 
to offer the image processing system which can perform distorted amendment at a 
high speed. 
[0009] 

[Means for Solving the Problem] (1) In the image processing system which creates the 
image after amendment from the image before amendment including a geometric 
distortion which the image processing system of this invention photoed the phantom 
which has specific geometric structure by the photography system, originated in the 
property of the photography system concerned, and was produced A measurement 
means to measure the coarse strain distribution which shows the distorted point of 
the specifying point specified according to the geometric structure of said phantom 
from the photography image of said phantom, A calculation means to compute the fine 
strain distribution which divides the image after said amendment into two or more 
minute fields, and shows the distorted point of the representation point in the minute 
field concerned using said coarse strain distribution, A concentration value decision 
means to determine the concentration value of the minute field on the image after 
said amendment based on the concentration value of the minute field on the image 
before said amendment specified by the distortion place of said representation point 
is provided. 

[0010] According to such a configuration, since the location of a pixel is associated 
for every minute field between the image before amendment, and the image after 
amendment in fine strain distribution, compared with the case where the location is 
associated for every pixel, the operation time can be shortened sharply. Therefore, 
distorted amendment processing is accelerable. This makes it the basis to consider 
that there is almost no distortion in the minute field. 



(2) The image processing system of this invention is equipment given in the above (1), 
and said calculation means is characterized by computing the fine strain distribution 
which shows the distorted point of the representation point concerned by carrying out 
the 2-dimensional interpolation operation of the coarse strain distribution which 
shows the distorted point of at least four specifying points near [ in said minute field ] 
the representation point. 

(3) The image processing system of this invention is equipment given in the above (1), 
and said calculation means is characterized by providing a means to divide the image 
after said amendment into two or more minute fields to which the core of an image 
differs in magnitude from a periphery, and a means to compute the fine strain 
distribution which shows the distorted point of the inner each representation point of 
said minute field based on said coarse strain distribution. 

A setting means for the image processing system of this invention to be equipment 
given in the above (1), and to set up an area of interest on the image after said 
amendment is provided further. (4) Said calculation means A means to divide the 
image after said amendment into two or more minute fields from which magnitude 
differs the inside of the field equivalent to the area of interest set up by said setting 
means, and out of the field, It is characterized by providing a means to compute the 
fine strain distribution which shows the distorted point of the inner each 
representation point of said minute field, based on said coarse strain distribution. 
The image processing system of this invention is equipment given in the above (1). (5) 
And said concentration value decision means A means to pinpoint the minute field on 
the image before said amendment equivalent to the minute field on the image after 
said amendment according to the distortion place of the representation point 
computed by said calculation means, It is characterized by providing a means to 
determine the concentration value in the minute field on the image after said 
amendment based on the concentration value in the minute field on the image before 
said amendment specified by this means. 

The image processing system of this invention is equipment given in the above (1), and 
said phantom has grid-like structure. (6) Said measurement means Distance during 
the lattice point of said phantom on the image after said amendment is integer-ized, 
and the coarse strain distribution which shows the distorted point of each lattice point 
when the integer-ization concerned was made is measured. Said calculation means It 
is characterized by dividing the image after said amendment so that it may be located 
at the core whose representation point of two or more of said minute fields is a pixel. 

(7) The image processing system of this invention is equipment given in the above (6), 
and is characterized by providing further a means to expand or reduce the image after 
said amendment by changing the distance during said lattice point. 

(8) The image processing system of this invention is equipment given in the above (7), 
and it is characterized by changing the distance between grids into the nearest 



integer. 
[0011] 

[Embodiment of the Invention] Hereafter, the operation gestalt of the image 
processing system of this invention is explained, referring to a drawing. Drawing 1 is 
the block diagram showing roughly the X-ray-diagnosis structure of a system 
concerning the operation gestalt of the image processing system of this invention. 
The system concerned is a system for presenting the diagnosis by a medical 
practitioner etc. with the image which performed the image processing which includes 
distorted amendment to the X-ray imaging image of the analyte obtained from the 
roentgenography system, and was obtained by this. 

[0012] It has the roentgenography system which consisted of X-ray TV cameras 104 
which picturize the image intensifier (I.I.) 103 which a diagnostic X-ray system detects 
X-ray source assemby 102 containing X-ray tube 100 and a beam limiting device 101, 
and the X-ray beam which exposure was carried out to Analyte M through the beam 
limiting device 101 from X-ray tube 100, and penetrated the analyte M concerned, and 
changes this into an optical image, and the optical image outputted from LI.103. and 
output the video signal according to this as shown in drawing 1 . The video signal 
outputted from TV camera 104 is stored in the image memory 206 which is changed 
into a digital video signal by A/D converter 201, and contains the image storing unit 
204 and the incidental information storing unit 205. 

[0013] Image data is stored in the image storing unit 204, and the incidental 
information which shows the photography conditions of whenever [ patient 
information /, such as a patient's identifier, sex, a birth date and weight, / , various 
information / accompanying the image data concerned /, i.e., distance between II-X 
line foci, dilation ratio, II size, camera station, and angle-of-coverage ] etc. is stored in 
the incidental information storing unit 205. 

[0014] An image processing system 202 performs various image processings including 
the distorted amendment by the distorted amendment section 203 to the X-ray image 
of the analyte read from the image memory 206. Distortion to which the image 
immediately after outputting from A/D converter 201 originates in the property of 
I.I.103 is included. The distorted amendment section 203 of an image processing 
system 202 processes the image data concerned according to the treatment process 
which reads the image data before amendment from the image storing unit 205 of an 
image memory 206, and is mentioned later, and creates the image after amendment 
with which the above-mentioned distortion was amended by this processing. 
[0015] The image data after the created distorted amendment is outputted to the 
image display section 207. The image display section 207 outputs the image data after 
distorted amendment to a monitor 210 or the laser imager 21 1. A display is presented 
with the X-ray-diagnosis image of the analyte after distorted amendment in this way. 
[0016] In addition, a computer (not shown [ Central Processing Unit : ]) manages 



actuation of A/D converter 201, an image processing system 202, and the image 
display section 207. Moreover, the image data of the image after distorted amendment 
can also be recorded on the magnetic disk (or optical disk) 209 which has large 
capacity with the instruction of a computer. In addition, both the incidental 
information stored in the incidental information storing unit 205 on the occasion of 
image recording is also recordable. 

[0017] In order to amend distortion which originated in the property of I.I. 103 and was 
produced, the distorted amendment section 203 of an image processing system 202 
performs distorted amendment processing as follows. That is, the distorted 
amendment section 203 is created according to the procedure which shows below the 
image after the amendment by which the distortion concerned was amended from the 
image before amendment including a geometric distortion which originated in the 
property of a roentgenography system and was produced. 
[0018] (1) Measurement of coarse strain distribution (measurement means) 

(2) Calculation of fine strain distribution (calculation means) 

(3) Decision of a concentration value (concentration value decision means) 

The measurement means of the above (1) measures the coarse strain distribution 
which shows the distorted point of the specifying point specified according to the 
geometric structure of a phantom on the image after amendment from the 
photography image of a phantom. The calculation means of the above (2) divides the 
image after amendment into two or more minute fields, and computes the fine strain 
distribution which shows the distorted point of the representation point in the minute 
field concerned based on the coarse strain distribution measured by (1). And the 
concentration value decision means of the above (3) determines the concentration 
value of the minute field on the image after amendment based on the concentration 
value of the minute field on the image before the amendment specified by the 
distortion place of a representation point. 

[0019] Thus, according to distorted amendment of this constituted operation gestalt, 
since the location of a pixel is associated for every minute field between the image 
before amendment, and the image after amendment in fine strain distribution, 
compared with the case where the location is associated for every pixel, the operation 
time can be shortened sharply. Therefore, distorted amendment processing is 
accelerable. This makes it the basis to consider that there is almost no distortion in 
the minute field. 

[0020] Here, the various parameters concerning both image before distorted 
amendment (henceforth a "distorted image") and image after amendment (henceforth 
an "amendment image") are defined concretely, and the distorted amendment 
processing concerned is explained to a detail. 

[0021] First of all, an image coordinate is defined as shown in drawing 2 . A coordinate 
(i, j) expresses the two-dimensional location of the pixel on the digital image inputted 



into an image processing system 202 (pixel). Let the unit of i and j be a [pixel]. For 
example, as image size shows in this drawing by the image of NxN, the pixel (1, N), and 
at the lower right of an image is expressed [ the pixel at the upper left of an image / 
the pixel (1,1), and at the upper right of an image ] for the pixel (N, 1 ), and at the lower 
left of an image as (N, N). At practical there is much what set image size to 
"1024x1024." Or there are some which are being set to "512x512" or "2048x2048." 
[0022] Next, the concentration value of a distorted image and an amendment image is 
defined as follows. 

D (i, j) — Concentration value C of the pixel in the coordinate (i, j) of a distorted image 
(i, j) — It supposes that a wire uses the grid arranged in the shape of a grid as the 
concentration value, next the above-mentioned phantom of the pixel in the coordinate 
(i, j) of an amendment image, and a grid coordinate is defined as shown in drawing 3 . 
(b) of a grid projection image in case (a) of drawing 3 does not have distortion, and 
drawing 3 shows the grid projection image in the case of having distortion. 
[0023] Next, the coarse strain distribution measured from a grid projection image is 
defined as follows. 

G~ (u, v) — Image coordinate on the distorted image of the lattice point (u, v) (two 
dimensions) 

the direction component of (Width j) of the image coordinate on the distorted image of 
the direction component Gj [ of (Width i) ] — (u, v) lattice point (u, v) of the image 
coordinate on the distorted image of Gi — (u, v) lattice point (u, v) — here, (u, v) are 
the grid coordinates showing the location of the lattice point, i.e., the intersection of a 
wire and a wire, here, it considers as the zero (iO and jO) of distortion amendment of 
the lattice point nearest to the core of a distorted image, and the zero of a grid 
coordinate is shown in (b) of drawing 3 according to this zero — as (uO and vO) — ** 
— it carries out. In addition, the above "G~" shall express a vector. 
[0024] For example, when the lattice point is set to 5x5, as shown in drawing 3 (b), (1, 
5), and the grid coordinate of the lower right lattice point are expressed [ (1 1), and 
the grid coordinate of the upper right lattice point ] for (5, 1), and the grid coordinate 
of the lower left lattice point as (5, 5), upper left coordinate, i.e., grid coordinate, of the 
lattice point. 

[0025] A certain lattice point (u, v) exists in a location which is different on a distorted 
image and an amendment image, respectively (image coordinates differ), and the 
correspondence relation of the image coordinate between both images expresses 
strain distribution. In other words, if distortion originally does not arise, the lattice 
point should be projected on the location of the lattice point on an amendment image 
(u, v), but when distortion mixes in I.I., a result projected on location G~ (u, v) of the 
lattice point on a distorted image (u, v) is brought. 

[0026] Next, the fine strain distribution computed based on coarse strain distribution 
is defined as follows. Correspondence relation between the predetermined sample 



point (representation point) on the amendment image computed from coarse strain 
distribution G~ (u, v) and the point (distorted point) equivalent to the representation 
point concerned on a distorted image is made into P~ (i, j), and is defined as follows. 
[0027] P~ (i, j) — Image coordinate of the distorted point of the sample point on an 
amendment image (i, j) (on a distorted image: two dimensions) 

Pi — (i, j) The direction component Pj of (Width i) of the image coordinate of the 
distorted point of the sample point on an amendment image (i, j) — (i,j) The sample 
point predetermined [ above-mentioned ] which is the direction component of (Length 
j) of the image coordinate of the distorted point of the sample point on an amendment 
image (i, j) is a sample point in the predetermined grid pattern on an amendment image. 
[0028] Measurement of the coarse strain distribution of the above (1) is specifically 
performed as follows. First, a grid is stuck on the front face of I.I.1 03, roentgenography 
is performed, and the image coordinate of all the lattice points (intersection of a wire 
and a wire) is searched for from the photography image obtained by this. Next, it is 
based on the dilation ratio of I.I.103, and spacing of the image number of dots for 
character on soft copy and a ****** wire, and asks for the spacing L of the wire on a 
photography image (the unit of a decimal is sufficient as [pixel [ ]]:L). 
[0029] If measurement of coarse strain distribution is restricted when there is no 
modification in whenever [ installation / of environments, such as earth magnetism, 
and a roentgenography system /, and angle-of-coverage ] etc., it is not necessary to 
measure it each time at the time of roentgenography (at the time of an image 
processing), and you may make it use the measurement data which were measured in 
advance and saved in the image memory 206. 

[0030] In addition, as a phantom which has geometric structure, it is not limited to the 
grid used here. For example, you may make it measure a small ball using the phantom 
arranged at equal intervals two-dimensional. About this measurement approach, the 
publication of bibliography "Beth Schueler et al. f "Correction of image intensifier 
distortion for 3D X-ray angiography", Proceedings of SPIE s Medical Imaging, and 
Vol.2432(Physics of Medical ImagingXpp 272-279 (1995)" is made to reference. 
[0031] Measurement of the fine strain distribution of the above (2) is specifically 
performed as follows. Here, based on the measured coarse distortion distribution, the 
pixel which shows the representation point in the minute field of the plurality in an 
amendment image computes whether it is equivalent to which pixel on a distorted 
image. 

[0032] The concentration value of C (i, j) is calculated in the following procedures. As 
shown on an amendment image first of all at drawing 4 (a), the grid grid pattern of LxL 
[a pixel] is considered from an image core (iO, jO). 

[0033] Next, in an amendment image, it samples at equal intervals mostly, and two or 
more sample points are acquired. At this time, a sample point is set as a pixel center 
position. Drawing 4 is drawing for explaining such sampling processing. Drawing 4 (a) 



shows the grid pattern (inside of a dotted line) chosen as the general drawing and 
arbitration of a grid grid pattern, and drawing 4 (b) shows the enlarged drawing of the 
selected grid pattern. In drawing 4 (b) f a black dot expresses a sample point, and the 
field shown by hatching expresses the minute field which makes the sample point 
concerned a representation point (representation pixel). By the way, since the 
distance during the lattice point is generally a decimal, the boundary of a grid does not 
pass along the core of a pixel. In such a case, it is made for the pixel which a bordering 
connoisseur requires to belong to the grid field of the larger one of the area divided by 
the boundary. However, since the form (magnitude) of a grid changes delicately with 
locations in such a case, it is necessary to change the method of division of a minute 
field for every grid, or to memorize the modification result in memory. These things 
cause the increment in the ******** processing time, and occupancy of memory. If 
the distance between grids is changed into the nearest integer to such a problem, all 
the forms of a grid can be made the same. Therefore, division becomes possible 
similarly about all grids only by memorizing the method of the division into the minute 
field of one grid. Compared with the case where can also lessen necessary memory 
space by this while being able to shorten the operation time, and the distance during 
the lattice point is set as a decimal, it is dominance. 

[0034] The inside of a grid pattern is divided here and it is made for each pixel to 
belong to the nearest sample point. The pixel surrounded by the dotted line is divided 
so that it may become the closest to the sample point which exists at the core than 
other sample points. In the example of drawing 4 , it will divide like a dotted line. The 
pixel in one grid surrounded by this dotted line is divided so that it may become the 
closest to the sample point which exists at that core than other sample points. 
[0035] The strain distribution of a sample point (i, j) is searched for by the following 
procedures 1 - the procedure 6. Procedure 1: Calculate the amount of gaps from the 
zero used for processing. For example, in the example of drawing 5 , the gap (HO, j-jO) 
by the image coordinate from the sample point (i, j) of arbitration to a core is 
converted into gap by the grid pattern. 
[0036] 

deltau= (int) [(HO) /L} 
deltav= (int) {(j-jO) /L} — (1) 

When the following numeric value is forward, a operator (int) omits below decimal point, 
changes it into an integer, in a negative case, revalues below decimal point, and 
changes it into an integer here. For example, after conversion is set to -2 when said 
numeric value is -1.2. Procedure 2: Ask for the strain distribution data (four points) 
used for processing. In the procedure 1 , the amount of gaps by the lattice point at the 
upper left of the grid pattern which contains (i, j) from an image core (delta u, delta v) 
was calculated. Therefore, the following formulas can express the grid coordinate of 
the upper left lattice point as (u1 and v1). The grid coordinate of the point of the upper 



right the lower left, and the lower right is similarly set to (u2, v1), (u1 and v2), and (u2 

and v2). 

[0037] 

u1 =u0+deltau v1 =v0+deltavu2 =u1+1 v2 =v1+1 ~ (2) 

Procedure 3: Ask for the distorted point of a point (i, j) based on surrounding strain 
distribution data. For that purpose, first as shown in drawing 6 (a), it asks for the 
internal ratio of the point in a grid pattern (i, j) in the image longitudinal direction and 
the lengthwise direction, and it is mi about them, respectively. : ni, mj : nj It carries out. 
[0038] 

m j = {j- (jO+deltau-L)} /L ni =1-mi — (3) 
mj = 0- (jO+deltav-L)} /L nj=1-mj ~ (4) 

Drawing 6 (b) shows the grid corresponding to the grid of drawing 6 (a), and the grid 
(grid) coordinate of Points P, Q, R, and S and the image coordinate of P, R, Q, S, T, U, 
and V become as follows, respectively. 
[0039] 

Account Grid (grid) coordinate of P (Gi (u1 and v1) and Gj (u1 and vD) 

The grid (grid) coordinate of Q (Gi (u2 and vD and Gj (u2 and vD) 

The grid (grid) coordinate of R (Gi (u1 and v2) and Gj (u1 and v2)) 

The grid (grid) coordinate of S (Gi (u2 and v2) and Gj (u2 and v2)) 

The image coordinate of P (Pi (i and j1) and Pj (i and jD) 

The image coordinate of R (Pi (i1 and j2) and Pj (i1 and j2)) 

The image coordinate of Q (Pi (i2 and jD and Pj (i2 and jD) 

The image coordinate of S (Pi 02 and j2) and Pj (i2 and j2)) 

The image coordinate of T (Pi (i and jD and Pj (i and jD) 

The image coordinate of U (Pi (i and j2) and Pj (i and j2)) 

The image coordinate of V (Pi (i,j) and Pj (i,j)) 

By applying 2-dimensional interpolation (primary interpolation function use : Bi-Linear 
interpolation) here, it asks for the distorted point of a point (i, j) as follows. As shown 
in drawing 7 (a) and (b), it computes independently about the direction of i, and the 
direction of j. 

[0040] [1] Presume G~ (u2 and vD to strain distribution G~ (u1 and v1) measured from 

the grid projection image, and P~ (i and j1). 

Pi CiJI )=mi Gi (u2 ,v1 )+ni Gi (u1 ,v1 ) 

Pj =(i and jD mi Gj+(u2 and v1) ni Gj — (u1 and vD (5) 

[2] Presume G~ (u2 and v2) to measured strain distribution G~ (u1 and v2) and P~ (i 

and j2) like the above [1]. 

Pi (ij2 )=mi Gi (u2 ,v2 )+ni Gi (u1 ,v2 ) 

Pj =(i and j2) mi Gj+(u2 and v2) ni Gj — (u1 and v2) (6) 

[3] Presume P~ (i and j2) to presumed strain distribution P~ (i and jD and strain 
distribution P~ (i, j) to ask finally. 



Pi Oj)=mj Pi (ij2 )+nj Pi (ij1 ) 

Pj =(i. j) mj Pj+(i and j2) nj Pj — (i and j1) (7) 

Procedure 4: Since Pi (i,j) of the coordinate (Pi (i,j) and Pj (i, j)) of the distorted point 
and Pj (i,j) are not integers in almost all cases, they search for the sample point of the 
nearest image from this point, ip = (int) (Pi+(i, j) 0.5), jp = (int) (Pj+(i, j) 0.5) — (8) 
Procedure 5: Assign the value on a distorted image to the value of an amendment 
image according to a degree type, assuming the field which has the same magnitude to 
the field (hatching section of drawing 8 (a)) divided on the amendment image that the 
coordinate (ip and jp) of the distorted point is set up on a distorted image at the core 
(hatching section of drawing 8 (b)), and there is no distortion in the interior. Thereby, 
the concentration value of an amendment image is determined. 
C(i+deltai, j+deltaj) =D (ip+deltai, jp+deltaj) 

(- Inside of deltaj<=L/[ L/2<=deltai, and ] 2 or a minute field) — (9) 

Procedure 6: Repeat a procedure 2 - a procedure 5 about all the minute fields that 

cover an amendment image top. 

[0041] The amendment image C (i, j) called for in the distorted amendment section 
203 concerned as mentioned above is saved to the image data storage unit 204. When 
the distorted amendment by this operation gestalt is compared with the 1st [ by 
different technique from the distorted amendment concerned ], and 2nd example of a 
comparison in respect of the amount of operations here, the result of an operation of 
the whole image is as follows. In addition, N and the number of pixels of a minute field 
are set to M for the number of pixels of the whole image. 

1st example integer-ized of comparison 4N which searches for an account and 
distortion distribution by primary approximation, and asks also for concentration by 
primary approximation, 17 Ns of addition, 10 Ns of subtraction, 20 Ns of multiplication, 
and 4Ns of divisions and distortion distribution are searched for by primary 
approximation. 2nd example integer-ized of comparison 4N which asks for 
concentration by zero-order approximation, 1 4 Ns of addition, 6 Ns of subtraction, If a 
distorted amendment integer-ized 4 N/M [ of 14 Ns of multiplication, and a 4Ns of 
divisions and this operation gestalt ], addition 14 N/M, subtraction 6 N/M, 
multiplication 14 N/M, and division 4 N/M minute field is made into 3x3 pixels Since 
the amount of operations which distorted amendment takes drops to 1/9 of the 
amount of operations which the example of a comparison of the above 2nd requires, if 
it thinks simply, it will become accelerable [ 9 times as many as this ]. 
[0042] Here, the various modifications of the above-mentioned operation gestalt are 
explained. 

(Modification 1) He is trying to sample the whole image similarly with the 
above-mentioned operation gestalt. However, there are few amounts of distortion as 
an image core, and they become larger, as it goes to the outside of an image from 
there. So, near an image core, it is few in a sampling, namely, the amount of operations 



can be further reduced by setting up a minute field more greatly. With the 
above-mentioned operation gestalt, the minute field was made into 3x3 pixels within 
all grid patterns, however — ****** there are few amounts of distortion centering on 
an image core about this in the field of every direction 1/2 of an image (area 1/4 of 
the whole image) — a minute field — 5x5 pixels — carrying out — except [ its ] — 
usually — passage 3x — it may be 3 pixels. If it does so, about 15% of the amount of 
operations can be reduced to the case of the above-mentioned operation gestalt, and, 
as a result, about 1.18 times as many improvement in the speed as this can be 
attained. 

[0043] (Modification 2) He is trying to sample the whole image similarly with the 
above-mentioned operation gestalt. However, it is near an image core that information 
important for a diagnosis exists, and it is only sufficient for the outside to display 
relation with the part (for it to be able to set up as an area of interest by an ROI 
marker etc.) currently observed. In such a case, the discontinuity (a thing like block 
distortion like DCT) of some structures seldom poses a problem, but just pinpoints a 
location from a great portion of structure. So, near the image circumference, it is few 
in a sampling, namely, the amount of operations can be further reduced by setting up a 
minute field more greatly. With the above-mentioned operation gestalt, the minute 
field was made into 3x3 pixels within all grid patterns, however — ****** there is little 
information important (area 3/4 of the whole image) outside the field of every 
direction 1/2 of an image centering on an image core about this — a minute field — 
5x5 pixels — carrying out — except [ its ] — usually — passage 3x — it may be 3 
pixels. If it does so, about 50% of the amount of operations can be reduced to the 
above-mentioned operation gestalt, and, as a result, twice [ about ] as many 
improvement in the speed as this can be attained. 

[0044] (Modification 3) With the above-mentioned operation gestalt, the minute area 
size of the image before amendment and the minute area size of the image after 
amendment are assumed to be the same. However, compared with the image core, the 
image circumference is expanded in fact by the effect of spool form distortion. Here, 
the dilation ratio based on images is set to 1, and the distance (for example, the die 
length of the wire set in the front face of a detector) projected on 1 pixel centering on 
an image presupposes that the dilation ratio of the location of the arbitration in an 
image is expressed by on what pixel it is projected in the location of arbitration. In this 
case, a dilation ratio is expressed as a function of the distance from an image core, for 
example, becomes like drawing 9 . An axis of abscissa is the distance from an image 
core here, an image core is shown 0% and the image circumference is shown 100%. 
Like 2x2 [pixel2] in a minute field, when comparatively small, since the difference in a 
dilation ratio seldom influences an amendment image, it can fully respond with said 
operation gestalt, but if a minute field is enlarged in order to attain improvement in the 
speed, the effect by the difference in a dilation ratio will become large, and the 



discontinuity on human body structure will arise on the boundary of a minute field. 
About this, it is solvable by constituting an image processing system like drawing 10 . 
The image processing system 300 shown in drawing 10 consists of the affine 
transformation section 301 , the image memory section 302, and the distorted 
amendment section 303. 

[0045] Here, the case (the maximum dilation ratio is set to 1.6 in this example) where 
I.I. has an expansion property like drawing 9 is considered. First, the distortion image 
read from the image memory 206 is inputted into an image processing system 300. 
The affine transformation section 201 changes reduction percentage for every [ to 
0.95 to 0.60 times ] 0.05, creates the contraction image of eight sheets, and stores 
these contraction image and the original image in the image memory section 302 
temporarily. 

[0046] Next, although the distorted amendment section 303 processes almost like the 
1st operation gestalt, it especially performs a setup of the minute field to a distortion 
image top as follows here. That is, the inverse number of the dilation ratio in the point 
sampling (representation) chooses the distortion image nearest to said reduction 
percentage, and minute field data are led from the image. This means leading the data 
of the almost same magnitude (magnitude on not an image but real space) as said 
minute field from the image before amendment reduced (expansion) so that the 
difference in the dilation ratio of each minute field under the effect of spool distortion 
may be canceled. 

[0047] (Modification 4) If spacing L of a wire is made smaller than an actual value in 
the above-mentioned modification 3, the probability for the pixel which cannot 
contribute to an amendment image with the property (property of the approach of 
searching for by zero-order approximation in case a value is calculated correctly) of 
the distorted correction method concerning this invention in spite of existing on a 
distortion image to arise will become high. Therefore, it is desirable to set up the 
spacing L of a wire more greatly than an actual value. 

[0048] For example, based on an expansion property as shown in drawing 9 , it 
considers expanding the spacing L of a wire by 1.6 times the actual value. That is, in 
the affine transformation section 301, a dilation ratio is changed for every 0.05 up to 
1.05 to 1.6 times, the expansion image of 12 sheets is created, the data of a minute 
field are led from the image expanded by 1.6 times about the image core, and the data 
of a minute field are led from the original image about the image circumference. It is 
made to lead the data of a minute field about the image field between an image core 
and the image circumference from the image with which the dilation ratio corresponds 
most in consideration of the both sides of expansion by the spool distortion of I.I., and 
image expansion. 

[0049] Expansion of the amendment image mentioned above is an example to the last, 
and according to the same processing as this, it can create the amendment image 



which has the dilation ratio of arbitration. Moreover, expansion of an image is not 
limited to the implementation approach by this modification, for example, may carry 
out expansion conversion of the amendment image directly. 

[0050] (Modification 5) In the operation gestalt of the above 1st, as drawing 4 (b) 
showed, the minute field which divided the inside of a grid grid at equal intervals was 
set up, but this minute field may be set up so that between grid grids may be straddled, 
as shown in drawing 1 1 . In this case, since the fine distortion distribution on a grid 
boundary can be searched for with primary interpolation, when not different [ the 
distortion place data of the lattice point can be used as it is as fine distortion 
distribution, and ] from the case where the total of fine distortion distribution data is 
the 1st operation gestalt, improvement in the speed can be attained rather than the 
1st operation gestalt. 

[0051] This invention cannot be limited to the operation gestalt mentioned above, but 
can deform variously, and can be carried out. For example, in the above-mentioned 
operation gestalt, although the diagnostic X-ray system with which the X-ray imaging 
image of analyte is obtained by the roentgenography system was explained, you may 
be the X-ray computerized-tomography scanning equipment with which the 
roentgenography system was equipped with other modality from which this differs, for 
example, I.I. 
[0052] 

[Effect of the Invention] As explained above, according to this invention, the image 
processing system which can perform distorted amendment at a high speed can be 
offered. According to such an image processing system, it becomes possible to make 
it work at clinical places, such as the usual test routine. Therefore, the precision of 
processing results, such as image analysis and measurement, can be improved, and a 
clinical place can be provided with the information for performing an exact diagnosis. 



DESCRIPTION OF DRAWINGS 



[Brief Description of the Drawings] 

[Drawing 1] The block diagram showing roughly the X-ray-diagnosis structure of a 
system concerning the operation gestalt of the image processing system of this 
invention. 

[Drawing 2] Drawing showing the definition of an image coordinate. 
[Drawing 3] Drawing showing a grid projection image. 
[Drawing 4] Drawing showing a grid pattern. 

[Drawing 5] Drawing showing the location within the strain distribution to be used. 
[Drawing 6] Drawing showing the correspondence relation of a grid. 



[Drawing 7] Drawing showing the operation of the distorted point. 

[Drawing 8] Drawing showing the correspondence relation between an amendment 

image and a distorted image. 

[Drawing 9] The graph which shows the expansion property concerning a modification 
3. 

[Drawing 10] The block diagram showing the configuration of the image processing 
system concerning a modification 3. 

[Drawing 1 1] Drawing showing the division approach of the minute field concerning a 
modification 5. 

[Drawing 12] Drawing showing the distorted mode concerning the conventional 
example. 

[Description of Notations] 



101 


— Beam limiting device 


102 


— X-ray source assemby 


103 


— Image intensifier (I.I.) 


104 


— TV camera 


201 


— A/D converter 


202 


— Image processing system 


203 


— Distorted amendment section 


204 


— Image storing unit 


205 


— Incidental information storing unit 


206 


— Image memory 


207 


— Image display section 


208 


— Computer 


209 


— Magnetic disk 


210 


— Monitor 


211 


— Laser imager 
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